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ABSTRACT 


A semi-empirical  formula  has  been  developed  to  fit  the  data  of 
Cavanagh  and  Rockman  on  thermal  radiation  from  2 oz  charges  of  high 
explosive  at  sea-level  and  at  pressures  of  12. 5,  1 .Q , and  6.5  mm  of 
Hg.  The  formula  is  based  on  the  point-source  model  of  Taylor  and  the 
mathematical  development  of  Bethe.  The  formula  indicates  that  the 
thermal  radiation  from  a high-explosive  charge  varies  inversely  as  the 
atmospheric  density  and  is  proportional  to  the  3.5  power  of  the  total 
explosion  energy^  which  is  taken  as  the  heat  of  combustion  if  the  ex- 
plosion occurs  in  a chamber  containing  more  than  sufficient  oxygen  to 
completely  burn  the  explosive;  the  total  explosion  energy  is  taken  as 
the  heat  of  detonation  if  there  is  insufficient  oxygen  to  burn  the 
explosive . 


INTRODUCTION 


mental  measurements  of  the  thermal  radiation  from  56.7  gm  (2  oz.)  charges 
of  high  explosive  (Pentolite)  detonated  in  a vacuum  chamber  at  pressures 
of  760  (sea-level)^  12. k,  7.0^  and  6.5  mm  of  Hg.  Figure  1 presents  their 
data  for  the  total  thermal  radiation  emitted  by  the  explosive  charges  as 
a function  of  the  initial  pressure  in  the  vacuum  chamber.  The  large  in- 
crease in  total  thermal  radiation  as  the  initial  chamber  pressure  is  re- 
duced from  760  to  72.5  mm  of  Hg  is  particularly  striking^  as  is  the  fact 
that  the  ratio  of  the  total  thermal  radiation  yields  is  very  nearly  equal 
to  the  reciprocal  of  the  ratio  of  the  initial  pressures.  As  the  initial 
chamber  pressure  is  reduced  still  further,  however,  to  7.0  mm  of  Hg,  the 
total  thermal  radiation  drops  by  a factor  roughly  two.  A clue  to  this 

behavior  is  provided  by  the  work  of  Filler—  who  has  shown  that  the  total 
heat  energy  produced  by  many  explosives  (including  Pentolite)  is  much 
larger  when  the  explosive  is  detonated  in  air  at  sea-level  initial  pres- 
sure than  when  the  explosive  is  detonated  in  a nitrogen  atmosphere.  The 
difference  in  energy  output  is  due  to  burning  of  the  explosive  in  the  oxy 
gen  of  the  air.  When  an  explosive  is  detonated  in  the  presence  of  a larg 
er  than  stoichiometric  amount  of  oxygen,  the  heat  energy  output  is  the 
'heat  of  combustion  "(H.  C.);  when  the  explosive  is  detonated  in  nitrogen 
(presumably  detonation  in  a vacuum  produces  the  same  result),  the  heat 
energy  output  is  the  ''heat  of  detonation"  (H.  D.)  . ^Pentolite  has  an 
H.  C.  of  2500  cal  gm  and  an  H.  D.  of  1100  cal  gm 

The  next  step  is  to  develop  a mathematical  model  of  an  explosion 
which  wil^^allow  the  ^otal  thermal  radiation  to  be  computed.  The  work 
of  Taylor—'  and  Bethe—  provides  the  basis  for  the  model. 


If  one  assumes  the  explosion  to  take  the  form  of  an  isothermal, 
expanding  spherical  shell  of  emissivity  one,  the  equation  for  the  total 


THEORY 


thermal  radiation  from  the  explosion  can  be  written  as 
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where  Y is  the  radius  of  the  spherical  shell 


a is  the  Stefan-Boltzmann  constant 


T is  the  temperature  of  the  spherical  shell 


Y is  the  initial  radius  of  the  shell,  and 
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Y is  the  radial  velocity  of  the  shell  (the  dot  over  the 
quantity  signifies  differentiation  with  respect  to  time) 


Now  the  total  kinetic  energy  the  expanding  spherical 


shell  is  obviously 
where  p is  the  density  of  the  undisturbed  air^  and 


(2) 


M is  the  mass  of  the  explosive, 
o 
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Bethe—  shows  that  E is  related  to  the  total  explosion  energy 

KJ.N 

Etot  equation 


^TOT  \lN 


(y  - 1)"^  + 0.807 


(3) 


y is  approximately  the  ratio  of  the  specific  heats  at  constant  pressure 
and  volume  of  the  surrounding  air  (y  is  1.4  for  air  at  room  tempera- 
ture and  normal  atmospheric  pressure) . 


Next  Y " can  be  expressed  in  terms  of  Y^  Y , M , ^ ^ and  y, 

by  use  of  Eqs.  (2)  and  (3),  as  o o 
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(Y  - 1)“^  + 0.807j  4TT3-^p(Y^  - Y^)  + 2~^ 


(4) 


To  express  T as  a function  of  Y,  it  is  assumed  that  T varies 
inversely  as  the  sum  of  and  the  mass  of  the  air^  swept  up  in  the 

expanding  spherical  shelly  that  is 


T T‘ 
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1 + 4TT3-^  p(Y^  - Y^)  i 

o o J 


(5) 


where  T^  is  the  initial  explosion  temperature. 

To  relate  T to  E„„„,  it  is  assumed  that  T is  proportional  to 

„ 4-  ■ ° TOT-^  o 

Etot^  that  IS 


T = C“^ 
o V TOT  o 


(6) 


is  a constant  having  the  dimensions  of  a specific  heat. 

If  the  values  of  Y~^ , T^  and  T , from  Eqs.  (4) , (5)^  and  (6)^ 
respectively^  are  substituted  in  Eq.  (1),  is  found  to  be  given,  for 

infinite  Y,  by 


RAD  o 

where  A = 0.5|^(y-1)"+  0.807 


= a C~^A°'^  0.4  (E^^^  p"^ 


TOT  o 


(7) 


-1 


2 


It  is  clear  from  Eq.  (7)  that  is  inversely  proportional  to 

the  initial  density  of  the  air  in  the  vacuum  chamber^  and  is  proportional 
to  the  3.5  power  of  the  explosion  energy. 

It  is  now  further  assumed  that  the  total  explosion  energy  is  given 
by  the  H.  C.  when  the  vacuum  chamber  contains  more  than  sufficient  oxy- 
gen to  completely  burn  the  explosive,  and  is  given  by  the  H.  D.  when  the 
amount  of  oxygen  in  the  chamber  is  much  less  than  stoichiometric. 


COMPARISON  OF  THEORY 
AND  EXPERIMENTAL  DATA 


2/  -1 
Filler—  |j.ves  the  H.  C.  of  Pentolite  as  2500  cal  gm  and  the  H.  D. 

as  1100  cal  gm  . Filler’s  vacuum  chamber  was  used  by  Cavanagh  and 
Rockman  for  their  experiments;  it  had  a volume  of  2050  liters.  At 
S.  T.  P.  this  volume  of  air  contains  613  gm  of  0^ • Take  the  charge 
weight  M as  an  even  57  gm;  then  if  it  is  assumed  that  the  combustion 
of  Pentolite  simply  burns  H to  H2O  and  C to  C0„  ) it  can  be 
shown  that  57  gm  of  Pentolite  requires  approximately  32  gm  of  atmos- 
pheric O2  for  complete  combustion.  Furthermore,  even  at  a vacuum 
chamber  pressure  of  72.5  mm  of  Hg,  the  chamber  still  contains  about  58 
gm  of  0„  which  is  more  than  the  stoichiometric  amount  (32  gm)  required 
for  complete  combustion.  At  a pressure  of  7.0  mm  of  Hg  in  the  vacuum 
chamber,  however,  the  amount  of  0„  in  the  chamber  is  only  5.6  gm,  much 
less  than  the  stoichiometric  quantity. 


It  seems  reasonable,  therefore,  to  use  the  H.  C.  of  Pentolite  in 
computing  the  total  explosion  energy  E for  chamber  pressures  of 
72.5  mm  of  Hg  or  higher,  while  the  H.  D.  of  pentolite  is  used  for 
chamber  pressures  of  7.0  mm  of  Hg  or  lower.  If  this  procedure  is  follow- 
ed in  computing  the  total  radiated  energy  E.  by  Eq.  (7)  and  if  the  values 
of  the  constants  A and  C in  Eq.  (7)  arecaken  to  make  the  computed 
^RAD  with  the  experimental  value  of  E^^  at  normal  atmospheric 

pressure,  then  Table  I and  Fig.  1 show  that^^^rae  computed  values  of 
E„._  at  pressures  of  72.5,  7.0,  and  6.5  mm  of  Hg  fall  within  about  157o 
orthe  experimental  values  of  Cavanagh  and  Rockman. 


Furthermore,  if  |he  ratio  o|  the  ^pecific  heats, 
1.4,  a as  0.562  - -----  — “ ^ - - 
as  1.292  X 10' 


y,  is  taken  as 

X 10  "2  ^fbt  cm  *"  , and  the  atmospheric  density 

gm  cm  — (0°C),  the  value  of  the  initial  explosion  shock 

front  temperature 


T can  be  computed  by  substituting  Eq.  (6)  into  Eq.(7) 
and  using  the  experimentally  measured  value  for  E^^^  as  given  in  Table  I. 
This  computed  value  of^T  is  about  8,100°K,  which  agrees  in  order  of 
magnitude  with  Erode 's—  computer  calculation  of  11,000°K  for  the  initial 
shock  front  temperature  for  a TNT  explosion  jn  air  at  s|a-level.  (TNT 
has  a much  higher  H.  C.  than  Pentolite,  3220—  cal  gm  . It  is  interest- 
ing that  the  ratio  of  11,000  to  8,100  is  roughly  equal  to  the  ratio  of 
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the  H.  C.  of  TNT  to  the  H-  C.  of  Pentolite.) 


CONCLUSIONS 


The  formula  derived  in  this  paper  gives  the  total  thermal  energy 
radiated  by  a high-explosive  charge  as  a function  of  the  total  explo- 
sion energy ; the  atmospheric  density^  and  other  quantities.  The  formula 
appears  to  be  useful  on  the  basis  of  the  very  limited  amount  of  data 
available.  Much  more  experimental  data  is  required  to  establish  the 
validity  of  the  formula  with  any  certainty. 

The  theoretical  justification  of  the  formula  is  tenuous  and  the 
formula  at  this  time  should  be  regarded  as  semi-empirical.  In  particular, 
the  mathematical  procedure  of  letting  the  shock  front  radius  Y become 
infinite  to  derive  Eq.  (7)  is  an  ad  hoc  stratagem  applied  to  fit  the 
experimental  data. 
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TABLE  I. 


Experimental  and  computed  values  of  the  total  thermal  energy  radiated 
by  high-explosive  charges  detonated  in  a vacuum  chamber  at  sea-level 
and  reduced  pressure. 


Pressure^, 
(mm  of  Hg) 

Thermal 

Yield^  Exptl.^ 
no  lens^ 
(cal) 

Thermal 

Yields  Exptl.^ 
with  lens^ 
(cal) 

Thermal 
Yields  9°^^” 
puted‘^ 
(cal) 

760 

4.49  X 10^ 

4.40  X 10^ 

4.40  X 10^ 

72.5 

5.23  X 10^ 

4.57  X 10^ 

4.61  X 10^ 

7.0 

2.48  X 10^ 

2.32  X 10^ 

2.70  X 10^ 

6.5 

3.01  X 10^ 

No  Data 

2.91  X 10^ 

a Cavanagh  and  Rockman  used  a lens  to  limit  the  field  of  view  of  their 
radiation  detector  in  some  of  their  measurements. 


b The  computed  data  are  normalized  to  the  experimental  value  of  the 
thermal  yield  with  a lens  at  760  mm  of  Hg. 
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FIGURES 


1.  Experimental  and  computed  values  of  the  total  thermal  radiation 
from  57  gm  Pentolite  charges  as  a function  of  pressure. 
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